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INTRODUCTION 


RATIONALE 

There will be relative motions between multiple Solar 
Power Satellites due to solar wind, lunar gravity, etc. 

This motion will cause beat frequency Interference on the 
earth's surface. Although the energy of such beats would be 
extremely small, the frequency of these beats might cause 
them to be of concern. 

BACKGROUND 

A Solar Power Satellite moving in a geostationary orbit 
appears to remain fixed with respect to any point on the 
earth's surface. However perturbations in that orbit due 
to solar wind, lunar gravity, etc., will cause the satellite 
to be in motion relative to a point of interest on the earth 
and the radial component of its velocity relative to that 
point will cause a Doppler shift in the frequency of the 
energy from the power beam as observed at that point. Assum- 
ing that the Reference retrodirective phase control system 
is being employed wherein the satellite on-board reference 
phase is derived from the pilot beam, there will be an 
additional Doppler shift. This is due to the radial com- 
ponent of the satellite's velocity with respect to the 
pilot beam source. A multiple number of satellites in motion 
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relative to each other and the earth will then produce energy 
at a multiple number of frequencies, resulting in a very 
large number of beat frequencies. 


SUMMARY OF RESULTS 


This report presents the results of an analysis of the 
frequency spectra and corresponding beat frequencies created 
by the power beams of a multiple number of Solar Power 
Satellites. 

METHODOLOGY 

Based on the assumptions given '<elow, the results were 
derived mathematically and verified through computer simu- 
lation. Frequency spectra plots were computer generateu. 
Detailed computations were made for the seven following 
locations in the continental USA: 


Houston, TX 
Seattle, NA 
Miami, FL 
Chicago, IL 
New York, NY 
Los Angeles, CA 
Barberton, OH 


GENERAL ASSUMPTIONS 


For the purpose of analysis, the frequency spectrum 
and beat frequencies created by Doppler shifts in the power 
beam frequencies were analyzed using the following general 
assumptions: 
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* 60 SoTar Power Satellites 

* Satellites located every 0.5° in longitude from 85.5° W 
(to 115°) 

* Each satellite is in an elliptically perturbed geo- 
stationary orbit as shown in Figure 1. 

— Each orbit has an eccentricity of 0.001 

— Each orbit has zero inclim tion with respect 

to the equatorial plane. (Note: The assump- 

tions concerning the shape and inclination of 
the orbits does not affect the general nature 
of the results of this report.) 

* The specific daily times that the satellites are at 

the perigees of their respective orbits ("perigee times") 
are statistically independent random variables. 

* Retrodirective phase control is used, except where 
otherwise noted. 

* The Rectenna (pilot beam source location) for each 
satellite is located within t 10° longitude of that 
of its associated satellite. 

* Equal sidelobe power is received at the seven geo- 
graphical locations. (This assumption was necessary 
because the power at those locations is primarily 

a function of the power in the far sidelobes of the 
power beams which are undetermined at this time. Jue 
to random phase errors, etc., at the space antennas, 
the far sidelobe patterns are li)cely to be determined 
primarily by the array factors of the individual power 
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FIGURE 1 


PATH OF SATELLITE IN ECCENTRIC ORBIT RELATIVE TO 
SATELLITE IN CIRCULAR SYNCHRONOUS ORBIT 
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modules whose broad beamwldths would cover the con- 
tinental USA. As soon as specific far sidelobe 
antenna patterns become available, their power levels 
can be incorporated into the results. 

FREQUENCY SPECTRA OF SATELLITE SIGNALS 







j 


The eccentric motion of a Solar Power Satellite causes a 
sinusoidal variation in the Doppler shift of signals trans- 
mitted to and from the satellite. The result is that the 
satellites *s power beam signal at the earth varies sinusoidally 
in frequency about the intended frequency of operation. For 
any of the 60 satellites and seven locations selected, we 
show that an eccentricity of 0.001 will result in a maximum 

freqnecy deviation, Af , of about 50 Hz in the power beam 

max 

frequency. This frequency deviation is proportional to the 
eccentricity. About half of the frequency deviation measured 
at any of the seven locations is due to the motion of the 
satellite with respect to its pilot beam source at the 
Rectenna and the remainder is due to the motion of the 
satellite with respect to the location studied. The period 
of the sinusoidal variation is 24 hours. It turns out that 
this result, which is shown in Figure 2, is analogous to 
what would result if the power beam from a satellite in geo- 
stationary orbit were frequency modulated with a sinusoidal 

signal whose "amplitude" and period were Af and 24 hours, 

Tnax 

respectively. 
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FIGURE 2 

POWER BEA.«1 FREQUENCY VARIATION DUE TO ECCENTRIC ORBIT 
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FREQUENCY IN HZ FROH 2.R5 GHZ 

FIGURE 3 

TYPICAL FREQUENCY SPECTRUM OF 60 SATELLITE POWER BEAMS AT HOUSTON - MIDNIGHT 
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FIGURE ^ 

TYPICAL FREQUENCY SPECTRUM OF 60 SATELLITE POWER BEAMS AT HOUSTON - 3 a.m. 



The 60 satellite power beams give rise to 60 signals 
sinusoidally varying in frequency, all with about the same 
maximum frequency deviation, but with random initial phases 
in the sinusoids of variation. A frequency spectrum that 
could exist in the Houston vicinity at midnight for a typical 
set of random perigee times is shown in Figure 3. The spectrum 
three hours later is shown in Figure 4. 

The distribution of power in contigous frequency bands, 
which we call the "Instantaneous Power Spectral Distribution 
(IPSD)", is shown in Figure 5 for the Houston vicinity at 
midnight for the perigee times used in Figures 3 and 4. The 
contiguous frequency bands are approximately 2 Hz wide. The 
IPSD is also equivalent to the percentage of frequencies that 
fall within each of the contiguous bands. The average of 
the IPSD's at midnight and 3 a.m. is shown in Figure 6. If 
we include in the average the IPSD's for 6 a.m. and 9 a.m. 
also. Figure 7 results. Including one IPSD for each hour, 
on the hour, starting at midnight gives the average shown 
in Figure 8. As more and more IPSD's are included, the 
average approaches the limit shown in Figure 9, which we 
call the "Average Power Spectral Distribution (ASPD)". An 
APSD can be obtained by continuous averaging of the IPSD's 
over a 24 hour period. It is symmetrical about 2.45 GHz. 

An APSD is equivalent to the spectral distribition that 
would result from 60 signals each frequency modulated by 
an appropriate amplitude sinewave for the case where each 



-2S c *n 

iM «i iM] 

nsyi! 7 

»mm V TK mwi^T, 5 *.*•.. « *.«.. mc i «.«. 
ilSTMnaitWS PMI VtCTMl DlST«IIUTlOi5 fO» HOUSIW 


-ft 0 *ft •» 

m MI fi«i 2. <5 M 

riu« a 

anutt V 2« 1RS^MT«(QUS Wl SncTUl simiurTiois, 

M fUF CACM NOUa Oi Tit «U>. fm HOUSTW 




niMMCt m Ml (MM ?.*S •>: 

MttK * 

uimi mca sKtraai aisTnuriOi w Nasrn 

(2 HI - HIK COVTlttOCS HttJOKy MSS' 


FIGURES 5-9 


OKIGINAL PAGE h 
OF POOR QUALm 


POWER SPECTRAL DISTRIBUTIONS - HOUSTON 



modulation index (ratio of maximum frequency deviation to 
modulation frequency) is very large. Figure 10 is the 
same APSD as Figure 9, except that greater frequency detail 
in the power distribution is shown by having the contiguous 
frequency bands be approximately 0.03 Hz wide (as compared 
to approximately 2 Hz in Figures 5 through 9) Figure 11 
is an expanded frequency scale version of Figure 10 to show 
the detail in the neighborhood of (2.45 x 10® + 50) Hz. APSO's 
for Seattle, Miami, Chicago, New York, Los Angeles, and 
Barberton, Ohio are shown in Figures 12 and 13. 

Note that the results shown in Figures 10 through 13 are 
nearly identical. At each location, the average distribution 
in frequency of the power due to the 60 satellites is horse** 
shoe-shaped and is concentrated in a band approximately 
100 Hz wide centered on the power beam frequency, 2.45 GHz. 

The reason that the APSO's are somewhat independent of the 
locations is that the radial velocity of a Solar Power 
Satellite with respect to any point in the continental USA 
is approximately the same. This is because the angle sub- 
tended by that portion of the earth's surface, as seen from 
each satellite, is small, due to the distance of the satellite 
from the earth compared with the earth's radius. 


•It is shown in the Detailed Analysis Section that the APSD 
is proportional to the power spectral density quantized over 
the contiguous frequency bands. 
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FIGURE 11 

AVERAGE POWER SPECTRAL DISTRIBUTION - HOUSTON 
(EXPANDED FREQUENCY SCALE, 2.A5 GHz + 50 Hz REGION) 
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AVERAGE POWER SPECTRAL DISTRIBUTIONS - VARIOUS LOCATIONS IN U.S 
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AVERAGE POWER SPECTRAL DISTRIBUTIONS - VARIOUS LOCATIONS IN U.S 
(EXPANDED FREQUENCY SCALE) 


BEAT FREQUENCIES 


The sum of two signals forms a composite signal whose 
maximum amplitude fluctuates ("beats") at a rate equal to 
the difference in the frequencies of the two signals (the 
"beat frequency") . This is shown in Figure 14 for the case 
of two signals of equal power, where f and f 2 are the 
Instantaneous frequencies of the two signals, and I fj^ ~ f£ I 
is the Instantaneous beat frequency. Thus power beam signals 
differing in frequency at the earth's surface due to differences 
in Doppler shifts will generate beat frequencies. 

Since the frequencies of the satellite power beams each 
vary sinusoidally with time, the beat frequency due to any 
two satellites varies with time over a frequency range from 
zero to some maximum frequency, with the period of the beat 
frequency variation being 12 hours. The maximum beat frequency 
can be from near zero to approximately 100 Hz depending 
primarily on the respective Initial phase of the sinusoids 
of variation in the two power beam frequencies. Figure 15 
illustrates two cases. For 60 satellites, there are almost 
1800 possible beat frequencies at any time. 

It is essential to recognize in Figure 14 that there 
is actually no signal at the beat frequency or at any other 
frequency except f^^ and f 2 » the instantaneous frequencies of 
the original two signals. The mere addition of two signals 
varying sinusoidally in amplitude at two different frequencies 
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AMplUiidt of Compotlto S<on«l. V|2 



• and V2 are two power beam signals of equal amplitude 
(power) with Instantaneous frequencies f^ and f2« 

• envelope of the maximum amplitude of the 
composite signal. 

• Uz the geometric mean of f^ and f2« 

• Ki ■ ^2! frequency. 

FIGURE 14 

BEATS BETWEEN TWO POWER BEAM SIGNALS RECEIVED >.'ITH EQUAL POWER 





cannot create a signal at a new frequency. Signals at the 
beat frequency can only be obtained through distortion of 
the waveform (of the two additive signals) by some non-linear 
process. One such process is thermal heating ^ which is 
proportional to the square of the instantaneous amplitude 
of the signal shown in Figure 14, and produces energy at 
the frequency Ifj^ - f 2 l» Any non-linearity which causes the 
generation of the product of two signals will produce a beat 
frequency signal. A non-linoar process which squares the 
composite (the sum) of the two signals does this — a "square- 
law" non-linearity. For the purposes cf this report , then, 
we assume that any non-linear process present is square-law . 

A square-law process is not only very reasonable to expect, 
but also one that allows us to conveniently calculate and 
look at the power spectrum of the beat frequency signals — 
the Instantaneous and the Average Power Sprecrral Distributions 
of the beat frequency signals (IPSDg and APSDg) . 

FREQUENCY SPECTRA OF BEAT FREQUENCY SIGNALS 

Consider a possible frequency spectrum for the 60 satellite 
power beams received, for example, at Barberton, Ohio at 
midnight. This is show-’ in Figure 16. For equal power 
received from all 60 satellites, and for 0.024 Hz-wide 
contiguous frequency bands. Figure 16 then represents the 
IPSD at midnight for Barberton. These 60 satellite signals 
generate almost 1800 beat frequency signals. Figure 17 
shows the frequency distribution of power for the beat frequency 
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signals in narrow contiguous frequency bands — the Instantaneous 
Power Spectral Distribution of the beat frequency signals 
(IPSDg) at midnight for Barberton, Ohio. The beat frequencies 
extend almost to 100 Hz. Figure 17 also represents the 
percentage of beat frequencies that fall within each of the 
contiguous bands. 

Since the time of day that each satellite is at its 
perigee is assumed to be a random variable, the power spectral 
distribution at any time can be expected to average out to 
be as shown in Figure 18 — the Average Power Spectral Distri- 
bution of the beat frequency signals (APSP^) for Barberton, Ohio 
The maximum frequency is approximately 100 Hz and the contiguous 
frequency bands are 0.03 Hz wide. Figure 19 shows that 
the APSD„ of Figure 18 is almost a linear function when 
plotted with a logarithmic frequency scale. 

The APSDg, as shown in Figure 18 and 19, is the statisti - 
cal average (statistical expectation) of the IPSDg, and is 
not necessarily the same as the 24-hour time average. How- 
ever, the time average of the IPSDg's should be very close 
to the APSD„ . 

D 

Figure 18 and 19, are for ali practical purposes, represent 
ative of the plots of an APSD^ at any location in the con- 
tinental U.S.A. 
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MITIGATING STRATEGIES 

Separate Transmit Frequencies For Each SPS 

The beat frequencies that have been discussed so far 
K are from retrodirective phase controlled satellites trans- 

mitting at nominally the same frequency and fall within the 
’ frequency range from zero to 100 Hz or so. By purposely 

separating the transmit frequencies of the power beams, 
the beat frequencies can all be made to fall above the zero 
to 100 Hz range because the minimum beat frequency will be 
the minimum frequency difference between any two power beams. 
However, to avoid any pumping in the ionosphere, that is, 
to avoid generating RFI due to the non-linear action of the 
ionosphere on microwave signals differing in frequency, the 
power beam frequency separation must be on the order of 
100 MHz. Therefore for a large number of power satellites 
we would run out of spectrum for the transmitting frequencies. 

All Transmit Frequencies Nominally the Same 

Since the beat frequencies for power satellites all 
transmitting at nominally the same frequency are caused by 
, differences in Doppler shifts, several mitigating strategies 

were investigated which minimize the differences or 
minimize the effects of the differences. These strategies 
fall into two categories — stationkeeping techniques and 
techniques that alter the transmit frequency of the power 
beat at the satellite. 
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stationkeeping - Maintain Satellites in Nearly Exact GEO 


f 


One obvious stationkeeping technique is to keep all 
satellites exactly or nearly exactly in GEO to essentially 
eliminate any relative motion of the satellites. Since a 
satellite's relative velocity is proportional to orbit 
eccentricity, then the Doppler shifts will be also. This 
means that since the maximum possible beat frequency is 
approximately 100 Hz (from Figure 18 and 19) for orbit 
eccentricities of 0.001, the orbit eccentricities would 
have to ben an order of magnitude better (less than 0.0001) 
to keep the maximum possible beat frequency to under 10 Hz. 
This technique is probably not practical and probably would 
require an exce«5Sive amount of fuel for stationkeeping. 


Stationkeeping - Synchronized Satellites 


A second stationkeeping technique allows for eccentric 
orbits, but requires that the orbit inclinations and eccentri- 
cities be identical and that the satellites be synchronized 
in such a manner that the magnitude and "phase" of the 
perturbations from GEO for each be identical at all times. 

This is equivalent to saying that at any time all satellites 
will be at the same relative point in orbits identically 
perturbed from GEO. Since the velocities of the satellites 
towards the center of the earth would therefore be identical 
at all times, the Doppler shifts at any point on earth due 
to the motion of the satellites will tend to be the same and 




' 
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will thereby tend to cancel out frequency differences* For 
the general assun^tions of this report, neunely zero orbit 
inclination, all eccentricities equal to 0.001, and satellites 
located over a longitudinal span of approximately 30° from 
85.5° to 115°, the maximum beat frequency on earth will be 
under 5 Hz. Wl^iether this stationkeeping technique is even 
feasible is another question. 

Transmit Frequency Altering Techniques 

Techniques that alter the transmit frequency of the power 
beam at the satellite do so in order to compensate for the 
frequency shift at earth caused by Doppler shift. 

Transmit Frequency Alteration - Determined by Frequency Shift 
at Rectenna 

One transmit frequency altering technique shifts the 
transmit frequency at the satellite by an amount equal to 
the frequency shift measured at the Rectenna, that is, by an 
amount equal to the difference from 2.45 GHz in the frequency 

it 

of the received power beam. If a retrodirective phase control 

_ 

Note that the 2.45 GHz reference source at each Rectenna, 
even if only used to determine the pilot beam frequency of 
satellites that are "perfectly" in GEO and hence have no 
relative motions, must be calibrated or referenced to some 
common frequency source so as to have a frequency difference 
less than 10 Hz. If the maximxjm compensated frequency 
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system is being employed, the satellite transmit frequency 
is shifted by shifting the pilot beam frequency by the 
measured frequency shift at the Rectenna. If a ground-based 
phase control system is being employed, such as Interfere- 

A 

meter-Based Phase Control , the transmit frequency at the sat- 
ellite can be changed in one of two ways as a result of the 
measurement of the frequency shift at the Rectenna. It can 
be changed via a command link or can be frequency-locked to 
a frequency which is synthesized from a pilot signal trans- 
mitted to the satellite at some suitable frequency, say 
X 2.45 GHz. In this case, the pilot signal frequency is 
compensated by the power beam frequency shift measured at 
the Rectenna. 

The power beam frequency shift at the Rectenna can be 
measured by direct comparison of the power beam frequency with 
a local oscillator at 2.45 GHz. The power beam frequency 
measurement/ correction loop only needs to operate very slowly, 
difference between satellites due to any mitigating strategy 
is, say 10 Hz, then the reference sources might have to be, 
say, within 1 Hzl 

*J. H. Ott and J. S. Rice "An Interferometer-Based Phase 
Control System”, Final Proc eedings of the Solar Power Satel- 
lite Program Review, Lincoln, Nebraska, April 22-25, 1980, 
pp. 289-292. (Available from National Technical Infor- 
mation Service.) 
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This allows long measurement Integration times for highly 
accurate determination of the frequency shift at the Rectenna. 
The slow operation of the loop is allowed by the low rate 
of change in the Doppler-caused frequency shift in the 
received power beam which is at most on the order of one 
Hertz per minute. 

For correction purposes, it may not even be necessary to 
directly measure the power beam frequency shift on earth 
but only to infer it. The power beam transmit frequency 
can be altered in a predictive fashion just by knowledge of 
the orbital characteristics. The power beam frequency shift 
at the Rectenna is a sinusoidal function of time, with the 
maximum magnitude of the shift being proporational to ec- 
centricity. The shift is zero when the satellite is at 
its apogee and perigee, and is at its positive and negative 
maximums 6 hours respectively from the time it is at its 
apogee and perigee. 

Transmit Frequency Alteration - By Means of 0n-3oard Atomic 
Clock 


Another transmit frequency altering technique uses an 
on-board atomic clock at each satellite to generate a refer- 
ence signal at a frequency of 2 x 2.45 GKz. When retrodirective 
phase control is being employed, the 4.9 GHz reference signal 
would be used as the reference signal for the phase con- 
jugation of the received pilot beam. Thus the phase con- 
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jugation process would automatically shift the power beam 
frequency by a compensating amount equal to the Doppler 
shift between the pilot beam source (the Rectenna) and the 
satellite. If ground based phase control is being employed, 
an uplink signal from the Rectenna would serve the same 
purpose as that served by the pilot beam in retrodirective 
phase control for generation of the power beam transmit 
frequency. 

All of the mentioned techniques that alter the transmit 
frequency of the power beam in a compensating manner result 
in a perfectly corrected frequency shift only at the Rectenna. 
The corrected frequency shift any place else is strictly a 
function of the geographical location of that place. Figure 
20 is a geographical plot of the contours of corrected 
frequency shift for a satellite located at longitude 100° W. 
These are maximum magnitude contours. For a satellite 
located at other than 100° W, the contours change longitude 
by an amount equal to the difference between 100° W and 
the longitude of the satellite. To find the maximum possible 
beat frequency between two satellites, assume that the 
two satellites are located in their respective orbits such 
that at a geographical point of interest, one satellite's 
corrected frequency shift is at its positive maximum when 
the other is at its negative maximum. The maximum possible 
beat, frequency is then just the sum of the frequencies of 









the two contours which exist at that point from the two 
satellites. From inspection of the contours of Figure 20, 
it can be determined that the maximum beat frequency that 
is possible (i.e., worst case) in the continental U.S.A. 
due to two satellites located anywhere from longitude 
85.5° W to 115° W is about 9 Hz. 





CONCLUSIONS 


Power satellites in orbits spanning 30*^ of longitude about 
longitude 100° W will generate beat frequencies from zero 
to a maximum frequency which is proportional to orbital 
eccentricity (all eccentricities equal) , and roughly pro- 
portional to the maximum eccentricity present (eccentricities 
unequal). In the continental U.S.A., for the above placed 
satellites: 

1. Eccentricities all equal to 0.001 can result in 
beat frequencies from 0 to 100 Hz for retrodirective 
phase control systems. Maximum beat frequency 

will be under 10 Hz for eccentricities less than 

0 . 0001 . 

2. Mitigating strategies have been determined which 
which maintain the maximum possible beat fre- 
quency to under 10 Hz for eccentricities up to 
0.001 (maximum possible beat frequency roughly 
proportional to maximum eccentricity present) . 
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DETAILED ANALYSIS 


DOPPLER SHIFT 

The Doppler shift in the power beam frequency of a Solar 
Power Satellite as observed at a point P^^ on the earth's 
surface, relative to the frequency transmitted by the satellite, 
is given by the expression 

(1) Af - - f i 
c 

where Af is the observed frequency shift, f is the frequency 
transmitted by the satellite, c is the velocity of light, and 
V is the radial velocity of the satellite with respect to P^ 
in a direction away from P^. Equation 1 also expresses the 
Doppler shift in the frequency of the pilot beam as received 
at the satellite, with P^ representing the location of the 
pilot beam source (center of Rectenna) , and f the transmitted 
frequency of the pilot beam. Both v and c remain defined as 
before. The v which corresponds to the P^ of concern is 
given by Equation 2. This is derived in Appendix A based 
on the geometry of the situation, with the symbols as defined 
in Figure A-1. 


PRECEDING PAGE BLANK NOT FUMED 


(r_ + Ax) + Ay^ + Az^ + r_^ 

G - E 

-2rg{[(rg + Ax) cos A6 + Ay sin A6]co« 0 Az sin 

For a zero inclinded elliptical orbit with eccentricity 
equal to 0.001, Az and the numerator term involving v^ are 
zero, and Ax and Ay* can )>(' neglected with a resulting maximum 
error on the order of only about 0.1 Hz in the Doppler shift. 
This is equivalent to ignoring the pertur)>ation in the 
position of the satellite from exact ideal GEO. Equation 2 
then simplifies to the following: 

r- 

v_(l - — cos 4> cos A9) - v_ ( — cos ♦ sin A6) 

* ^ “G 

V ■ 


Figure 21 depicts the situation. From Appendix A, numerical 
expressions for v and v are given by: 

K 1 

(4a) Vjj ■ 3.07 sin ~ t (neters/sec. ) 

(4b) v^ « 6.13 cos ^ t (meters/sec.) 

*Figure 1 shows that Ax and Ay have maximum values of 42 km 
and 84 km, respectively. 
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where t is the time, in hours, of the satellite from its 
perigee, and T > 24 hours. Vj^ and are proportional to 
the eccentricity. 

It was assumed that a Rectenna site (pilot beam trans<* 
mitter location) in within + 10° longitude of its corres- 
ponding power satellite. Thus |A6| < 10°, and the Doppler 
shift in the received pilot beam frequency at the satellite 
can be approximated by letting v ■ in Equation 1. This 
is equivalent to 4* * 0, and A6 > 0 in Equation 3, which 
for calculation purposes, puts the pilot beam transmitter 
on the equator at the same longitude as the satellite. A 
more accurate expression for this v would be obtained by 
continuing to let A6 ■ 0 in Equation 3 and letting 4 equal 
some mean latitude of all the Rectennas. Ho ever, the 
maximum error in the Doppler shift using v ■ Vj^ was deter- 
mined by calculation to be only on the order of about 0.1 Hz, 
and therefore Vj^ is used in this report for the Doppler 
shift in the pilot beam fzequency. 

Let fp be the transmitted frequency of the v^lot beam, 
fj^ be that of the power beam at the satellite, and Vj^ and 
Vg^ be the velocities of the satellite w' th respect to its 
associated Rectenna (pilot beam transmitter), and with 
respect to P^, a location of interest on earth. The 
frequency shift, Af^, observed at relative to f^, due 
to Doppler shifts, is then given by 



Pq in geostationary orbit. 


FIGURE 21 

GEOMETRICAL RELATIONSHIPS AMONG THE EARTH AND 
SATELLITES IN GEOSTATIONARY AND ECCENTRIC ORBITS 
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DEFINITIOJJS OF SYMBOLS 


■ arbitrary point on the earth's surface 

P^ ■ instantaneous position of satellite relative to P^ 

P * mean position of satellite relative to Pp, equal to 
G t 

position of a satellite in geostationary orbit 
relative to 

^ - latitude of P£ is positive in northern hemisphere) 

A0 » longitudinal difference between P^, and P^. (A6 is 

positive if Pg is west of Pj.) 

Vg ■ instantaneous satellite velocity relative to P^ 

Vg «/v^ + v^ corresponding 

orthogonal components of Vg with directions defined 
as follows: 

Vj^ : radially away from center of earth 

v^ : direction of flight of geostationary satellite at Pq 

v„ : normal to equatorial plane (toward North Star) 

rj. * radius of earth 
r^ « radius of geostationary orbit 

d * distance between Pg and P^ (positional variation of Pg 
with respect to P^) 

d ■ /(Ax)^+ (Ay)^ + (Az)^, where Ax> Ay, Az are corresponding 
orthogonal components of d with directions respectively 
the same as those defining Vj^, v^, v^. 


FIGURE 21 (Continued) 

GEOMETRICAL RELATIONSHIPS AMONG THE EARTH AND 
SATELLITES IN GEOSTATIONARY AND ECCENTRIC ORBITS 


. * ■« 
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T 


(5a) tf - - V, 


SR 


C- ^SE 


(5b) 


f f 

V - -2. V 

C SR C SE 


(5c) 



+ V 


SE 


) 


where Vg^ is expressed by Equation 3. Substituting the 
expressions for Vj^ and from Equation 4 into Equations 3 
and 5 gives a numerical expression for which is proportional 
to the eccentricity and is of the form; 


(6a) Af^ = - A sin (|r- 1 - a) 


where A is the maximum magnitude of sind o is the numerical 

value of the phase of Af^ with respect to that of v^^ . 

Let us denote A as Af . Equation 6a can then be rewritten as; 

max 

(6b) Af = - Af sin (t - t ) 
e max T o 


T 

where t^ * ^ Irrespective of the location of in the 
continental U.S.A. , Af has a value in the neighborhood of 
50 Hz. The magnitude of t^, which is the difference in time 
between when the radial velocity reaches its maximum and when 
Af reaches its maximum, is less than half an hour. (The 
instantaneous frequency, f ^ , at is given by 


f = f + Af 
e o e 


f - 
o 


Af 


max 


sin 


2it 


(t - 


‘o>- 


This means that the power beam signal at can be represented 
by a sinusoidally frequency modulated signal 


. / 
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V » V sing (t) 


where 

B(t) = |(2irf^)dt, 

with a inoudlation period of 24 hovirs. Performing the indicated 
integration, with f^ obtained from Equation 7, yields 

V = V sinr2nf t + Af T cos (t - t ) + B 

L o max T o oj 

where 6^^ is an arbitrary phase. (If f^ and are in Hertz, 

then t and T are in seconds.) The product Af T is the 

in 3 X 

modulation index, the ratio of the maximum frequency deviation, 

Af , to the modulating frequency, 1/T.) 
max 

If t* is the time of day (0 _< t' £ 24 hours) , then 

(7) t = f - t^ 

where t^ is the time of day that the satellite is at its 

perigee. Combining Equations 7 and 6b provides the following 

form for Af ; 

e 

(8) Af = - Af sin |t' - (t + t, ) I • 

Since the modulation index is very large, there is power at 
the following discrete frequencies at a given time of day, 
t' , from 60 satellites: 

(9) f = f - Af sin ft* - (t + t, ) I ; n =1,2,... 60. 

en o max T L on In J 

n 

The t, are assumed to be statistically independent continuous 
In 

random variables, each with a uniform probability density 
function given by 


- 37 - 


(10) p(t^„) - ij, 0 i < 24. 

In other words, any has equal probability of being any 

time during the day. Figure 2 is a plot of Equation 9 for 

a particular n, with t' - (t + t. ) as the horizontal 

on xn 

axis variable. 

INSTANTANEOUS POWER SPECTRAL DISTRIBUTION 

For equal power received from all satellites, the power 
in a frequency band at a given time of day is proportional to 
the number of spectral lines in that frequency band. The 
distribution of power in all bands forms a histogram, which 
when expressed on a percentage power basis, defines the 
Instantaneous Power Spectral Distribution, IPSD. This 
was depicted in Figure 5. Performing a time average of the 
IPSD's yields the Average Power Spectral Distribution (APSD) . 
(It is sufficient to average over a period of 24 hours.) 

MATHEMATICAL EXPRESSION FOR THE AVERAGE POWER SPECTRAL 
DISTRIBUTION 

The APSD, expressed as a percentage distribution, is 
equivalently given by the expression 

60 [^1 

(11) APSD(f^) = ^ I S^(f)df, 

n-1 J 

^ 1 = all integers from 1 to + ® 

where S (f) is the power spectral density of the signal 
n 

received from satellite #n, P is the average power received 
from each of the 60 satellites, and f^ defines the frequency 
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bands, each band being Af wide, where 


(12) « Uf. 


The power spectral density of the signal received from 
a satellite with instantaneous frequency given by Equation 9 
is, from Appendix B, 


^ i ; I f - f I < Af 

^ /Af ° 

/ max o 


(13) S (f) 
n 


0 ; f - f > Af 

' o' max 


where f^ is the transmitted frequency of the pilot beam, and 
is the average received power. Using Equation 13, the power 
in a band Af wide at a frequency f^^ due the signal from 
satellite #n, with average power = P, is: 


( 14 ) 


f^+ Af 

S (f)df 
n 


° + 1 " 


- Af , or f . - f > Af 

max 1 o max 

n n 


£ Fsin'^ + "] , f £ < - 

7T Af 2 i o — 

L max J 

n 


A£ , and 
max 

n 


f, ^ - f >-Af 

i + 1 o — max 




sin-1 U- * 1' - sin-1 


Af 


max 


Af 


P 

IT 


7T 1 f. - f 

JL i O 

2 Af 


max 


f, - f I 

1 o I 
/ < Af 

max_ I — max 

n 

and 

1*1*1- *o 

i I'*... 

t, . , - f„ > if 5 

. 1 + 1 o — max , and 
' n 

f, • f < Af 
1 o — max 

n 


where f 


i + 1 


+ Af. 


Given Af and the Af , Equation 11 for the APSD can then be 

m A X 

n 

numerically evaluated using Equations 12 and 14. Typical 
APSD *8 for two different Af's were shown in Figures 9 and 10. 


MATHEMATICAL RELATIONSHIPS AMONG THE INSTANTANEOUS AND AVERAGE 


POV^ER SPECTRAL DISTRIBUTIONS AND THE POWER SPECTRA^. DENSITY. 


From Appendix C» 


, ou 

IPSD(f ) - i I 7(f 
^ oU ^ 




i = all integers from 1 to + “ 


where 


V(f, - f ) = 
i en 


1 at f . for f . < f < f + Af 
I i 1 — en i 

0, otherwise 


and where is defined by Equation 9, and f^ defines the 
frequency bands, each band being Af wide, where f^^ is des- 
cribed by Equation 12. Thus, from Appendix C, it is proven 


that 


APSD(f 


l’ ■ f 


IPSD(f^)dt 


is equivalent to Equation 11. Appendix C also develops the 


fact that 


APSD(f^) « E^j^UPSD(f^) } 


where E^j^{IPSD(f } is the statistical expectation with 
respect to the times in the day when the satellites can be 





at the perigees of their respective orbits. 


Let S(f) be the power spectral density due to the 60 satellites 
at a point of observation. Then, 


60 

(16) S(f) - I S„{f) 

n-1 " 


The average 
f^ is given 

(17) 

In terms of 
can then be 


of S(f) over a frequency interval, Af, at a frequency 


by 


S(f^) 


_ 1 _ 

Af 


f^i 

S(f)df, 


i = all integers from 1 to + ® 
Equation 16 and 17, the APSD given by Equation 11 
expressed at a frequency, f ^ , as: 


(18) APSD(f^) = 11^ ST?7) . 

As Af approaches zero, sTf^) approaches S(fj^) in value, and 
thus 

APSD(f ) 

(19) Lim — — i- » S(f, )/60P 

Af-^0 ^ 


where f^^ represents any frequency, and S(fj^)/60P is the 
power spectral density, normalized in terms of total power. With 
the property given by Equation 19, the APSD, as expressed 
by Equation 18, can be thought of as a power-normalized, 
quantized version of the power spectral density. 
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BEAT FREQUENCIES 


Two received satellite signals, (t) and V 2 (t) , of 
equal amplitude (power) at slightly different instantaneous 
frequencies, fj^ and f 2 » give rise to a resulting composite 
signal of the form: 


( 20 ) 

where 


v^2^t) * Ej^ 2 (t) COs27Tf^2t 


f 1 - f2 

E^2<t) » 2Vcos2tt ( ^ )t. 


" ^2 


12 


given that 

where 

(21a) 


V 12 (t) = Vj^ (t) + Vj (t) , 


Vj (t) = V cos 27lf^ t 


and 


(21b) V2(t) = V COS 


The composite signal, v^ 2 (t) (Equation 20), and the amplitude 
function, Ej 2 (t), are shown in Figure 14. The frequency, ^^ 2 ' 
is the geometric mean of f ^ and f 2 . 


Beats, the maximum in amplitude of v^^ 2 » occur at a rate 
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since the maxitnums occur every half period of E ^2 shown 
In Figure 14. We will only consider frequencies to be 
psoitive in value and since f^ ** negative, we will 

define the beat frequency to be 



Beat frequency signals are formed by the product of two 
signals. For example, for Vj^ (t) and V 2 (t) described by 
Equation 21: 

2 

Vj^(t) • V 2 (t) “ V cos 2 tt f^t cos 2irf2t 

= — cos 2tt (f^ - ^2^ ^ 

+ Y ~ cos 27t (f^ + f^) t 
where | - fj) 1 is the beat frequency. 

For a square-law non-linear process, 

[vj^(t) + V 2 (t)J^ = v^{t) + 2 Vj^(t)V 2 (t) + V 2 <t), 
and the beat frequency comes from the Vj^(t)V 2 (t) term. 

BEAT FREQUENCY EXPRESSION - SINUSOIDALLY VARYING FREQUENCIES 

At a given time of day, t* , power received at a point of 
interest on earth from the 60 satellites is at the 60 
sinusoidally varying frequencies given by Equation 9, and 
the beat frequency due to any two satellites which, is given 
by the magnitude of the difference in their frequencies, is: 
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( 22 ) 


iAf 


max 


^ .in [f - * t,,)] 


- Af sin ~ ft* - (t , + t, ,)] 

maxj T L oJ IJ J 

where i and j refer to any two satellites, the "i th" and j th" 
ones. Equation 22 yields an expression of the following form 
for the beat frequency: 

(23) |sin f - Yjj) I 

where is a phase term, and is the maximum value of 
the beat frequency. can have a value from 

when (t . + t, , ) * (t . + t, .), i.e., when the variation in 

01 11 Oj ij 

f , and f , are "in phase", to a value of lAf + Af I 
ei ej ' maxj max^ ' 

when 

<'01 * ‘li> - '‘oJ * ‘ij’ * 

i.e., when the variation in f^^^ and f^j are "180° out of 
phase." This was alluded to in Figure 15. 


Af 


nax. 


- Af 


max 


DETERMINATION OF TOTAL NUMBER OF BEAT FREQUENCY SIGNALS 


Sixty satellite signals yield 

60 ! _ 60(60 - 1)_ ,7..^ 

21 ( 60 - 2 ) 1 * 2 

possible beat frequency signals at any time. The satellite 
signal frequencies vary sinusoidally with time (Equation 9) 


44 


and the phases of these sinusoids n ■ 1,2, . , »60) 

are random. Therefore, since the period of a beat frequency 
variation is 12 hours, there will be an infinite number of 
possible beat frequencies over a range from zero to approxi>> 
mately 100 Hz during a 12 hour period. 

INSTANTANEOUS AND AVERAGE POWER SPECTRAL DISTRIBUTION OF 
BEAT FREQUENCY SIGNALS 


For equal power received from all satellites, the iPSDg 
is (from Appendix D) expressed by 


60-1 60 


(24) 


IPSDg(f^) = 


60 (60 - 1) 


£ I J'h - - * 

n»l n*m+l 


en 


), 


where 


i « ail integers from 1 to + * 


V(f. - 


f - f 1) 

em en ' 


1 at f. for f. < |f - f 1 < f + Af 

i 1 — ' em en ' 1 

0, otherwise 


and where f and f are defined by Equation 9, and f. defines 
the frequency bands, each band being Lf wide, where f^ is 
described by Equation 12. 


Appendix D shows that Equation 24 is equivalent to the 
IPSDg obtained via self-correlation of the IPSD: 


(25) IPSD, (f ) 
0 1 . 


2 * 60 
60-1 



j 


), 


i ■ all integers from 1 to + *. 
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i 


In general. Equation 25 is computationally ainpler than 
Equation 24. 

The APSD^ is defined as the statistical expectation of the 
IPSD^ with respect to the times in the day when the satellites 
can be at the perigees of their respective orbits: 


(26) APSDj(f^) » E^^{lPSD^(f^) } 


i » all integers from 1 to + ® 


Appendix D shows that Equation 26 is equivalent to obtaining 

the APSD_ via self-correlation of the APSD: 

B 


(27) APSD_(f.) -= I APSD(f jAPSD(f, - f J , 

® ^ ^ j.O ^ ^ ^ 

i » all integers from 1 to + ®. 
Also, from Appendix 0, analagous to Equation 11, 


rfi + Af 

(28) APSDg(f^) - Sg(f)df 


where S^(f) is the power spectral density of the beat frequency 
signals and can be obtained from the self-correlation of the 
power spectral density, S(f): 


(29) Sj(f) ' S(f)S(f - f)df 


where the power spectral density was defined in Equation 16, 


-46- 


MITIGATING STRATEGIES 


1. Altering the Power Beaiw Transmit Frequency to Compensate 
for Frequency Shift at the Rectenna 

Techniques that correct for the frequency shift at the Rectenna 
resulting from Doppler shifts by altering the transmit frequency 
of the power beam at the satellite all have the same corrected 
frequency shift at a point on earth. This is true whether 
direct or inferential ground measurement or phase conjugation 
with an on-board atomic clock is being usod. The maximum 
magnitude of the corrected frequency shift is a function of the 
geographical location of that point. The governing equations 
are developed below. 

a . Power Beam Frequency Shift Being Determined by Direct 
or Inferential Measurement . 

To the expression for the uncompensated frequency 
shift of a power beam at a point, P,, denoted as Af , 
we add a correction term which is the amount that the 
power beam transmit frequency is altered and is equal to 
the negative of the Af^ determined at the Rectenna. This 
yields the net or residual frequency shift in the power 
beam frequency at point P^ which we call the corrected 
frequency shift. Denoting the Af^ determined at the 
Rectenna, as Af>, and At* the corrected frequency shift 
at point Pg, we have 


-47 


(30) Af = Af - Af_. 

e e R 


1. ) Retrodirective Phase Control Being Employed 


The uncompensated frequency shift, is given 

by Equation 5c: 
f- 


Af 


- — ‘''r * ''SE> 


Therefore, since Af_ is that Af for the Rectenna, 

Re 


“r ■ ■ -F •''r ^ ''sr’ ' 


where is the Vg^ pertaining to the Rectenna. Using 
the previously established fact that Vgj^ = Vj^, 

(31) Af^ = - ^ v^. 


Substitution of the expressions for Af^ and Af^^ from 
Equations 5c and 31 into Equation 30 yields 


(32) |4f;| = ^ |v, - Vjjl. 


2.) Ground-Based Phase Control Being Employed , 


The uncompensated frequency shift, Af^, is given by 


(33) Af^ = - Vgg. 

Therefore Af , the Af for the Rectenna, is given by; 


(34) Af = - -2. V = - -2- v„. 

R C SR C R 


Substitution of the expressions for Af^ and Afj^ 
from Equations 33 and 34 into Equation 30 yields 
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which is identical to Equation 32. 

b. Power Beam Frequency Shift Being Determined by Phase 
Conjugation with An On-Board Atomic Clock 


Let be the Doppler shift at the satellite of a 
signal transmitted from the Rectenna, referenced to the 
frequency# the desired frequency of the power beam at the 
Rectenna. At the satellite, the pilot beam frequency (retro- 
directive phase control) or the synthesized pilot frequency 
(ground based phase control) is then given by; 


(35) f = f + Lf , 
o o o 

As a result the phase conjugation of the signal at frequency, 
f^, with a conjugating signal at frequency 2f^, generated by 
an atomic clocJc, the power beam transmit frequency, is 


(36) f = 2f - f = f - Af 
1 o o o o 


Let Afj^ be the Doppler shift in the power beam frequency at 

a point P_, on earth. Therefore the power beam frequency 

£ 


at Pg is 


(37) f f , + Af = f - Af + Af, 

el 1 o o 1 

Now Af^ is defined by the expression 


Af 

e 



f . 

o 


Therefore 

(38) |4f; I - |£^ - £^| > |i£j - AfJ. 
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Now 


(39) = - -2. v_ 

o C SR 




and 

(40) Af, 


V 

C SE 


C ''sE* 


Substitution of the expressions for Af^ and Afj^ in Equations 

39 and 40 into Equation 38 gives a result 
f 




identical to Equation 32 again. 


c. Evaluation of the Maximxim Magnitude of the Corrected 
Frequency Shift at a Point on Earth . 

The corrected frequency shift magnitude, |Af^|, expressed 

by Equation 32, has a maximum value, lAf I which can be 

' e ' max 

evaluated by using Equation 3 to express v__, where v_ and 

^ K 

v^ are given by Equation 4, and by using the fact that the 
maximum value of (A sin H' + B cos 4') equals /a* + B* , where 
A and B are arbitrary magnitudes, and V is an arbitrary phase. 
The result, in Hertz, is: 
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yield the contours of Figure 20 when the power satellite 
is located above the equator at longitude 100® W. (The 
signal horizon in Figure 20, obtained from numerical 
solutions to Equation A-24 in Appendix A, assumes that 
the power beam signals do not propagate over the optical 
horizon of the earth as seen from a GEO satellite at 
longitude 100° VJ.) Generation of contours for power 
satellites located at other than longitude 100° W can 
be used in conjunction with Figure 20 to show that the 
maximum beat frequency possible in the continental U.S. 
for orbit eccentricities equal to 0.001 is about 9 Hz. 

2 . All Satellites Synchronized - Perturbations From GEO 
Identical At All Times 

Two satellites which are in zero inclined orbits with 
eccentricities both equal to 0.001 and synchronized such 
that the magnitude and "phase" of their perturbations from 
GEO are identical at all times have identical relative 
positions at all times in orbits described by Figure 1. 
Therefore, for the two satellites, the velocities v^^ and 
v^ given by Equation 4, which we denote here as 
Vj ^2 ^Tl' respectively identical at all 

times, that is, 

(42a) Vjj E E Vg = 3.07 sin ^ t 

(42b) Vjj E E - 6.13 cos ^ t 
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Referring to Equation 5c, if and Af ^2 respective 

Af 's at a point, P-, on earth due to the two satellites, 

6 Cft 

then we can establish the relationships 


tf 


el 


C ^'^Rl '^SEl^ 


Af 


e2 


■ — '^SE2> ' 


where and are the two Vg^'s at point corres- 
ponding to the two satellites. Let fgj^2 beat 

frequency due to the two satellite power beams. Therefore, 


<«> *B12 = “el - ^ll- 


Since Vjj 5 and i 


*B12 c ''SEI - ''SE2 


★ 

Equation 44 can be solved for the maximum value of f^^^ ^^12^ 

using Equation 3 to express the Vg^^'s, using the relationships 

given by Equation 42, and remembering that the maximum value 

/ 2 ? 

of (A sin '1' + B cos 4*) equals / A + B , where A and B are 
arbitrary magnitudes, and S' is an arbitrary phase. The result, 

in Hertz, is: 


*The maximum beat frequency is identical to in Equation 
23 (for i = 1, j « 2) . 
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(1 - “ cos if' cos ^8.) 

G Z 

1 + (— ) - 2 — cos 4) cos A6, 

^G ^ 


-.*5 


( 45 ) 


(1 - — cos 4) cos A 02 ) 
G 


/l + (— ) * - 2 ^ 


cos 4» cos A6, 


Pl2 = 3.07 ^ 


r r. 


+ 2 ‘ 


— cos 4> sin LQ. 
^G 


/ 1 + (— ) - 2 - 


cos 4 cos A6, 
^G ^ 


— COS 4) sin A9^ 
r _ 2 


/ r 2 r 

1 + (— ) - 2 — cos 4) cos A6_ 

r r_ 2 

G G J 

where = 2.45 x 10^ Hz, c = 3 x 10® m/sec, and where A6j^ 
and A02 are the A6's corresponding to the two satellites. Now 
here, is a function of 4>» A0j^, and A02. Searching out the 
values of 4>, A0p and ^ that maximize F ^2 yields 

(46) Fj^2 ^ 5 


Therefore, the maximum beat frequency is less than 5 Hz when 
all the satellites are in orbits with eccentricities equal to 
0.001 and synchronized such that their perturbations from GEO 
are always identical. 


ORIGINAL PAGE Ib 
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appendix a 


APPENDIX A 


I. VELOCITY OF A SATELLITE IN AN ORBIT PERTURBED FROM GEO 


RELATIVE TO A POINT ON EARTH 


a. Radial VelocitN 


Referring to Figure A-1# a satellite, whose orbit is per- 
turbed from a geostationary orbit has a radial velocity, v, 
with respect to a point on earth, determined as follows. Let 

the rectangular coordinate system x , y , 2 be such that the 

000 

x^ - y^ plane is the equatorial plane, the North Pole lies on 

the + z axis, and the Prime Meridian lies in the x - z plane, 
o 00 

6g is then the longitude (East) of P^, and 0 = 90° - is its 

latitude ((p- is the colatitude). The radius of the earth is r„ , 

If i , j , k are unit vectors in the positive x , y , 2 
o -'o o 0 0 0 

directions, respectively, then the location of P^. is defined 
by the vector, r^, where 


(A-l) rj =- 


Xj = sin»j cos 


(A-2) i Yg = sin<|)p sin 6, 


- ■'e *£• 
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(Symbols defined on 
next page.) 


Flight path of 
satellite in 
geostationary orbit. 



FIGURE A-1 


GEOMETRICAL RELATIONSHIPS BETWEEN THE EARTH AND 
A SATELLITE IN AN ORBIT PERTURBED FROM A GEOSTATIONARY ORBT' 








DEFl:nTIONS OF SYMBOLS 


■ arbitrary point on the earth's surface 
P^ • instantaneous position of satellite relative to P^ 

P^ ■ mean position of satellite relativ.*^ to P^, equal to position 
of a satellite in geostationary orbit relative to Pg 
• latitude of P^ (positive in northern hemisphere, 
negative in southern) 

" « colatitude of P^ 

0^ ■ longitude of (positive in eastern hemisphere, 
negative in western) 

« longitude of (positive in eastern hemisphere, 
negative in western) 

V- ■ instantaneous satellite velocity relative to P (and to P ) 

Vg • / ^ v^, where Vj^, v^, are corresponding 

orthogonal components of Vg with directions defined 
as follows: 

Vj^ : radially away from center of earth 

: direction of flight path of geostationary satellite 

at P^ 

: normal to equatorial plane (toward North Star) 

V = radial velocity of satellite with respect to 
r^. * radius of earth 

r_ * radius of geostationary orbit 
tg = vector from center of earth to Pg 
r^g * vector from p£ to 

a « angular difference in direction between satellite velocity 
and radial component of satellite velocity with respect to P^ 
d » distance between Pg and P^ (positional variation of Pg 
with respect to P^ ) 

d * /TTx) ^ ^ (^y) ^ t i^z)\ where ix, ^y, Lz are corres- 
ponding orthogonal components of d with directions 
respectively the same as those defining v , v , v„ . 

R T N 

FIGURE A-1 (Continued) 

GEOMETRICAL RELATIONSHIPS BETWEEN THE EARTH AND 
A SATELLITE IN AN ORBIT PERTURBED FROM A GEOSTATIONARY ORBIT 
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iiet Pg be the instantaneous position of the satellite relative 
to Pg, with instantaneous velocity, v^. If Pq, the mean 
position of the satellite relative to P^, is the position of 
a satellite in geostationary orbit relative to P^, then the 
vector, d, which defines the deviation of the satellite from 
its mean location, represents the perturbation of P^ from 
geostationary orbit. For convenience, we define d in the 
rectangular coordinate system x^, yj^, centered on Pg. 
is in the direction of - r_, the negative of the vector from 
P- to the center of the earth, and y, lies in the earth's 
equatorial plane. The magnitude of r^ is the radius of the 
geostationary orbit, d is made up of orthogonal component 
vectors Ax, Ay, Az such that 

cT * Ax + Ay + Az 

If ij^, jj^, are unit vectors respectively in the positive 
Xi» y^^, directions, then 

(A-3) r^ = - ir^ 
and 

(A-4) ct » i^Ax + j^Ay + kj^Az 

The radial velocity, v, is the component of Vg in the direction 
of r__, the vector between P,. and P_, and its magnitude is 
expressed by 

(A-5) V ■ Vgcos a. 
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Equation A-5 is obtained by taking the vector **dot" product, 


(A-6) V - Vg • e^g, 


where e^g is a unit vector in the direction of r^g. 


fore, 


(A- 7) e. 


ES 


ESI 


Substituting Equation A-7 into A-6 yields 


(A- 8) V 


V * T 

S ES 


I 

I *7 } 


If r is a vector from the center of the earth to P, 

5 • 


(A-9) d » r^ + rg. 


Now 


(A-10) rg = (rg) = Vg 


Therefore, 


which becomes, using Equations A-3 and A-4, 


Vg - - i^r^ + i^Ax + jj^Ay + k^Ay. 


But 


r » 0 
G 


Therefore, 


There- 


, then 
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(A-11) Vj ■ jj^Ay + k^Az. 


(A-12a) ^ ijAx 


(A-12b) - j^Ay 


(A-12C) Vj, - kjAz 


Then Equation (A-11) becomes 


(A-13) Vj - ijV, jjVj + kjV^ 


(Note from Figure A-1 that Vj^ is the velocity component normal 
to the equatorial plane# and as Ax and Ay approach zero, relative 
to r_, v_ and v_ respectively become the components of v- 
radially away from the center of the earth and in the direction 
of the flight path of a satellite in geostationary orbit.) 


'es ' “"s ■ ’'e 


Using Equations A-1 through A-4 and Equation A«-9, Equation 
A- 14 becomes 




(A-15) 


- ®E 


+ sin e, 


+ k^r^ cos 0^) 
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Now it it teen that the unit vectors of the y^, and 

O' ^0 0 

Xi» coordinate syttemt are related aa follows: 


(A- 16) J 


•1 • io ®C ^ ®( 


®G ^ ®®* ®G 


)c - k 
1 o 


Subttituting the values for ij^, jj^, from Equation A-16 
into Equations A-13 and A-15, performing the operations 


indicated in Equation A- 3, and then letting 


90 - 4> and 

O 


(A-17) A6 - 6^ * 


results in the following expression for v: 


(A-18) v » 


Vjj[(rg + Ax) - rp cos ♦ cos A 6 j 


+ v^[Ay - r^ cos <{» sin A 6 ] 
+ Vjj[Az - rg sin 


(tg + Ax)* + Ay* Az* + r^* 


- 2rgi[(r^ + Ax) cos A 6 + Ay sin A 6 ] cos (> 

+ Az sin 41 } 


b. Relative Orbit Velocities v^ and v_ for Zero Inclination 

^ T 

Orbit with Ecce ntrici ty » 0.001 


The ellipse of Figure 1, as depicted in Figure 21, with 
parameters given by Figure A-1, is described by the equation 
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(A-19) 


2 


1 


Ax 


Ax 


o' + 


max 


'Ay 


max 


where Af and are the maximum magnitudes of Ax and Ay, 

max max ^ 

By substitution into Eqation A-19, it is easy to show that 


(A-20a) :x = - Ax cos ^ t 

ma X 1 


(A-20b) Ay * Ay sin —?■ t 
•' ■'max T 

provide a '•olution to Equation A-19, where T = 24 hours, and 
t is the time, in hours*.. r\ the satellite from its perigee. 
From Equation A-12, 



Therefore, ^aking the time derivative of Equation A-20 yields 

(A-21a) ''r = “f" T ^ 

(A-21b) ^ cos ^ t. 

The numerical values of Ax and Ay , which are found from 

max max 

inspection of Figure 1, are 

Ax = 42.164 km 

max 

Ay = 84.328 km. 

-'max 

Note from Figure 1 that Ax and Ay are proportional to the 
orbital eccentricity. Thus Vj^ and v,j are proportional to 
orbital eccentricity. Since T = 24 hours = 24 x 3600 sec.. 
Equation A-21 then becomes 
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(A-22a) 


Vj^ ■ 3.07 sin t (meters/sec. ) 

MV 

(A-22b) » 6.13 cos t (meters/sec.) 

II. OPTICAL HORIZON OF THE EARTH AS SEEN FROM A SATELLITE IN 
GEO 

The optical horizon of the earth, as seen from a satellite 
in GEO, establishes the limit on the earth's surface over 
which a power beam is assumed not to propagate. An example 
is that shown in Figure 20. This horizon is the locus of 
points at which the vector r^g in Figure A-1 is 90° to the 
vector r^ (with Ax * Ay * Az ® 0) . In vector "dot" product 
notation, this relationship between rj.g and r^ can be expressed 
as 


• ^ES = 0 

where 


■'e • ■'es = ■'e "^es ’ 


and y = 90 is the angle between r^ and r^g. 


Using equations A-1, A-2, A-15, A16, and A-17, Equation 
A-23 becomes (with Ax = Ay = Az = 0 in Equation A-15) : 

(A-24) CCS 4i cos 6 = — , 

G 

which yields the locus of points making up the horizon in 
terms of and A9 . 
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APPENDIX B 


POWER SPECTRAL DENSITY OF A RECEIVED 
SATELLITE POWER BEAM SIGNAL 


APPENDIX B 


POWER SPECTRAL DENSITY OF A RECEIVED SATELLITE POWER BEAM SIGNAL 


Consider a power beam signal whose frequency varies accord 

ing to Equation 9 which is repeated here for satellite In, 

with t' - (t + t, ) replaced by t": 
on In n 


(B-1) f * f - Af sin ^ t" 
' ' en o max T n 

n 


(Since t* is the time of day, t" varies from 0 to 24 during 
a T = 24 hour period.) Figure 2 is a plot of Equations 12 
and D-1. 


Now consider the energy received over a 24 hour period 
in a frequency band Af Hz wide as shown in Figure B-1. If P^ 
is the average power in the received satellite signal, then, 
since the frequency is varying so slowly, all of that power 
is essentially received in the frequency band from f to 
f + Af during the time intervals from t„ to t^ + At and 
from tj^ to t^ + At. Thus the energy in the band at frequency 
f^ is represented mathematically as 

(B-2) E^ (f^) = 2(P„At) 

But the energy in the frequency band is the power in the 
frequency spectrum from f to f + Af integrated over a 
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FREQUENCY 



POWER BEAM FREQUENCY VARIATION DUE TO ECCENTRIC ORBIT 
-"DWELL TIME" IN A FREQUENCY BAND 


24 hour period 


®.“.> 



+ Li 

S^(f) df dt 
n 



+ Li 

S^(f) df 


where T » 24 hours. Therefore, 


ff + Af 

T * S (f)df » 2 P Zit 
Jf n n 


or 



+ Af - 

S (f)df » ^ P 
n T n 


At 

-n 


2 p Ai 
T n Af 


As Af 0, 


(B-3) 



+ Af 

S (f)df 
n 


s (f ) 

n a 


Now t” can be considered to be a function of f so that 
n a 


Af 


t;; (f^ + Af) - t" (f ) 
n a n a 

n 


Therefore, as Af-^0, 

df 

a 

Since the frequency changes so slowly with time, Equation B-2 
is still valid because even as the change in frequency, Af, 
gets very small (Af 0), the corresponding time interval. At, 
is still very large. Therefore we can equate Equations B-3 
and B-4: 


(B-4) f P„ 


At 

- 1 P 

Af 

T 
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Since is an arbitrary frequency, 
(B-5) s„(£) -|p„ 


From Equation B-1: 


- h 


Now f varies in frequency, f, from f - f to f + Af 

an o max o max 

n n 

Therefore, 


T - 1 ^0 ” ^ 

<®-«' - 5 ? ft — ■ 

max 


f - f 


< Af 
— max 


Performing the indicated derivative operation in Equation B-5 
using the value of t^ and conditions on f from Equation B-6 
yields the power spectral density: 


(B-7) S„(f) 

n 


1 

TT 


/Afl - (f - f 

na X 0 

n 


f - f < Af 

o — max 


0; f - f > Af 
' o max 


APPENDIX C 

MATHEI-IATICAL DEVELOPMENT OF THE INSTANTANEOUS AND 
AVERAGE POWER SPECTRAL DISTRIBUTION IN TERMS OF POWER 
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MA THEMATICAL DEVELOPMENT OF THE INSTANTANEOUS AND AVERAGE 
POWER SPECTRAL DISTRIBUTIONS IN TERMS OF POWER SPECTRAL DENSITY 


From Equation 9, the instantaneous frequency of the "n th" 

satellite's power beam signal as received at a point on earth is 

f = f - Af sin ~ ft' - (t + t, )1 
en 0 T I on In J 

n 

where T' is the time of day. Consider a period of time, about 
a specific time of day, t^, so that 

(C-1) t' - t' << T = 24 hours. 

O 

Therefore the corresponding satellite power beam signal, , as 
received at that point on earth is 

(C-2) V = V cos [2Ttf (f - t') + B 1 
n n I en o nj 

where 


= maximum signal amplitude 


en 


- Af 

sin 

il 

t' - 

(t + 

max 

n 


T 

o 

on 

- Af 

max 

n 

sin 

2ir 

T 

t' - 

(t + 

on 


B^ = arbitrary phase • 


PRECEDING PAGE BLANK NOT FILMED 


- 75 - 


All of the & , (n « 1,2,3...N), (N « 60 tatellites) , are 
n 

considered to be uncorrelated (specifically > statistically 
independent) because the N satellite power beams are not 
phase locked together. 

Now a sinusoid at a frequency f^, with maximum amplitude, V, 
has its power spectral density, S, given by 



2 

where V /2 is the average power, and 5(f - f^) is an impulse 

A 

at f » f^. Therefore the power spectral density, S(f), of N 
satellite power beam signals measured over a short period 
time at a specific time, t^, during the day is given by 

^ N N 

(C-3) S(f) ‘I S„(f) * I P 

n-1 " n-1 " 

2 

where is the average power, V^/2, of a signal. (If two 

A 

or more are at the same frequency, we will consider S(f) 
to be the statistical expectation of the power spectral 
density at a time, t', with respect to the phase, This 

means in this case that the powers are additive at a 
frequency because the B^ are uncorrelated. However, the 
percentage of time during a 24 hour period when two of the 
60 sinusoidally varying frequencies are the same is, for all 
practical purposes, essentially zero) . 

We define the Instantaneous Power Spectral Distribution (IPSD) 
as 
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(C-4) IPSD(f^) - ^ 


f. + Lf 

X 


S(f)df, 


i ■ all integers from 1 + • 

where n is the normalization factor so that 


(C-5) I IPSD(f.) - 1 (100%) 

1-1 ^ 

where detines the frequency bands, each band being hf 
wide, where f^ is described by Equation 12. Substituting 

/s 

Equation C-4 into C-5, with S(f) in C-4 given by Equation C 
solving for n yields 


00 

[f, + Lf 

■ N 

n - I • 
i«lj 

1 

y p 6(f - f ) 

n en 

n-1 

» 4 


which becomes 


I p 


n«l 


6 (f - f ) df - I P 
en ‘'-I n 

n»l 


For this report, we assume that all P^^ are the same power, P 


Therefore 


and 


n « Np 


(C-6) IPSD(fj) - i J 


I £ (£ - £„) «. 

1 


3 and 
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Since 


+ tt 


6(f - f') df - 
n 


‘l i *.n ‘ *1 ♦ 
(O; otherwise. 


the IPSD can also be expressed as 


(C-7) 

where 


IPSD(£i) - i ! 7(f^ - f., 
" n«l 


(C-3) A(f 


^ *” (O; otherwise . 


The twenty- four hour time average of the IPSD is the Average 
Power Spectral Density (APSD) , and is given by the expression: 


(C-9) APSD(f^) - i 


IPSD(fj)dt'; T « 24 hours 


Substituting for IPSD(fj^) from Equation C-6 yields: 

>T 


APSD(fj) - i 


1 

H* 


'i V' 


T 6 (£ - f )df df 

*'i *n 

n«l 


(C-10) APSD(f^) 


NT 


N 

ff, + 


4 

rT 

I 

1 


M£ - f.„)df 

n -1 J 



0 


df 


Solving for t' from Equation 9, we see that t* can be thought 
of as a function of f , that is, 
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f (f ) - T? Bin'^ — + t + t- 

• n 2tt Lf on In 

nax 


Making a change In variables from time to frequency in the 
inner integral of Equation C-IO and taking the derivative 
of t’ (fgjj) with respect to f^^^ then yields 




f (0) *" ^ 

en 


df 


^f^ - f ) 


«n 


nax 


en o 


where 


and 


f (0) « f when t* » 0 
en en 


f (T) « f when t' ■ T. 
en en 


But from Equation 13 we see that 


1 

IT 


n en 


A 


- ‘o> 
n 


Therefore Equation C-^10 becomes 

APSD(fi) - hi I ^ J, 

n * X n X « 


(C-ll) 


‘.n“» 


Now 


79 - 


(C-12) 


f (t) 

•n 

6 (£ 


• n 


■ 2S„(f), 

n cn •n n 


because, as can be seen from Figure B-1, f 


•n 


f , twice , for 


f (0) < £ (f ) < £ (T). 
ti — n n 


Substituting the results of Equation C-12 into C-11, and 
remembering that ■ P gives 


(C-13) APSD(f^) 


1 ^ 

HP 


N 


I 

n*l 


S„(f)df 


But Equation C-13 is the same as Equation 11 (N ■ 60) . There- 
fore, if the APSD is thought to be defined by Equation C-9, 
then we have just proved that Equation 11 is a valid equivalent 
definition. Note that Equation C-13 is equivalent to 


(C-14) APSD(f^) 


where 


HP 


I S(f)df 


(C-15) S(f) 


N 

n»l 


Note that imbedded within the inner integral in Equation C-10 
is the fact that 


(C-16a) S (f) 

a 



S (f)dt' 
n 
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(C-16b) S(f) 




because, from Equation C-3, with * P, 


S„;f) = P6 (f - f ). 
n en 


Now the statistical expectation of S(f) with respect to the 


t, is given by 


Eti {S(f)} = [ 

^ C 


• • S (f ) p ^^2^2 1 o • • 


IN' — 11 "12 


Since the t^j, are statistically independent. 


p(tii,ti2, . . .,tijj) - p(t, j^)p(t^2^ • • 

Using this fact and, from iJquation 10, the fact that 


(C-17a) Ej.^{S(f)} = I S (f)} 

n*l ” 


(C-17b) 


= Ie {s„(fD 

n»l 


(C-17c) 


= I i (f)dt, 

n-1 T Jo “ 


(C-17d) 


1 rT /V 

--4 s 





because 


(C-18) i^(f) . PM£ - 

that iSf S (f) is a function of only one of the t, not of 
any other. From the form of Equation 9, we seet that we 
can interchange t* and t^^^ in the derivation subsequent 
to Equation C-10 to enable us to equate the appropriate 
expressions in Equation C-16 and C-17 to yield 

(C-19a) E^, {S (f) } = (f) 

tin n 

and 

(C-19b) E^, {S (f) } = (f) 

tin n 

Specifically, Equation C-19 says that the statistical expect- 
ation of the power spectral density at any time, t', with 
respect to the time of day, ®tc, that the satellites 

can be at their perigees, is the sane as a 24 hour time 
average of that power spectral density. This is tantamount 
to saying that the statistical expectation of the IPSD equals 


the APSD. 
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SPECTRAL DISTRIBUTION OF BEAT FREQUENCY SIGNALS 

Beat frequency signals are formed by the product of two 
signals. The product of two signals results from square-law 
processing of the linear sum of two signals. Consider the 
square-law process which then has the input-output relationship 


out 


= KvJ . 
in 


If V symbolizes v ^ when 
omn ■' out 


then 


V. * V + V , 
in m n 


V = X(v + V ) ^. 
omn m n 


When and are described by Equation C-2, with 

t" = t' - t << T = 24 hours, 
o 


(v + V ) ^ ® + 2v V 

m n m n m n 


(D-1) 


2V* + 2V^ cos 
n n 


(V 2 V ' 
(2 *2 


cos 


V 2 V 2 

+ -^ + -B. cos 
2 2 


27T(2f 

)t" 

+ 2B 

en 

n 

2iT(2f 

)t" 

+ 26 

en 

m 


+ 2a 


= n 


+ (continued on next page) 
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(D-1 

+ V V cos 
m n 

’2TT(f + 

em 


t" + 

m 

»n>] 


cont' d. ) 

+ V V cos 
m n 

'2TT(f 

em 


t" + 

<®m “ 
ID 

®n>] 

|j . m m. 


Only the last tern in Equation D-1 is in the audio frequency 

range and is generally not at D.C. (zero frequency). It is 

generated by the product This is the beat frequency 

signal due to two signals and is denoted here by where 

v„ is that portion of v in the frequency region | f - f | , 
Bmn omn em en 

that is. 


(D-2) 


Bmn 


= KV 


m n [ 


2tt (f 


em 




(f - 




+ (6 


m 


- s„> 


)' 


m / n. 


Therefore the power spectral density, Sg^^(f), of this beat 
frequency signal measured over a short period of time at a 
specific time, t^^, during the day is given by 

(KV V ) ^ 

S„ (f) = — 6(f - (f - f._)); f > 0, m ^ n. 

Bran 2 em en 


Let So (f ) be the totality of the power spectral density of 

beat frequency signals from all N = 60 satellites around 

t' = t*. Therefore, 
o 

(D-3) Sj(f) - I I ■” 

m«l n»l 


Now 


f 

en 


f 

em 


Therefore 


(f. 


em 


f 

en 


) 
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^ (KV V )* 

* —■ F- « ' f > 0, m ^ n. 


Suppose Then 6(f + (f^„ - f.„)) only has a 

en en cm ©n 

value atf“-(f -f )<0. But f > 0. Therefore 

em en 

✓V ^ 

S„ (f) « 0, and S- (f) is non-zero. Therefore, in 

Bnin DiDn 

/V ^ 

general, either or will be zero, and Equation 

D-3 becomes 


(D-4) 


SB(f) 


N-1 N 

I I 

m»l n*m + 


V ) 

ID n 


6(f 



f 

em 


) 


(If there are more than one beat frequency signals at a fre- 
quency, for example, |f^^ - f^^ | = |f^y - then we 

/N 

will consider Sg(f) to be the statistical expectation of the 

power spectral density at time, t', with respect to the phase 

6 . This means that the cowers are additive at that fre- 
n ■ ■ ■■ — 

quency because, as discussed in Appendix C, the are 
uncorrelated.) 


We define the Instantaneous Power Spectral Distribution 
of the beat frequency signals, IPSD„, as 

D 


(D-5) IPSDg(f^) = - 


f, + Af 

X ^ 


Sg(f)df; 


i = all integers from 1 to + ®, 


where n is a normalization factor so that 
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(D-6) I IPSD.(f.) « 1 (100%) 


Substituting Equation D-5 into D-6, with Sj(f) in Equation 
C-5 described by Equation D-4, and solving for n yields 


n 


2K* 


N-1 V. 

I i p p 

m«l n-m+1 " ” 


where 

(D-7a) 

and 

(D-7b) 


m 



V -2 



For this report, we assume that all are the same power, P. 

Therefore, P and P both equal P. Therefore, 

in n 


n = 2(KP) -ilNzl) „ (KP) * N(N - 1) 


and 


(D-8) IPSDj(fj) - 


'l ♦ ' 

N-1 N 

I I «(f - |£ 

m»l n«o+l 


en 


f |)df 
en ' 


Equation D-8 can also be expressed 


(D-9) 


IPSDg(f^) 


2 

N(N-l) 


N-1 

I 


m*l 


N 


1 V(f 

n*m+l 



f 

en 


) 


where the V function is defined in Equation C-8, 
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Now with IPSD(fj^) given by Equation C-7, consider the 


self-correlation of IPSD(f^) 


A * 

IPSD(f. )*IPSD(f J * I IPSD(fJ IPSD(f, - f.) 

^ j«0 ^ J ^ 


’ Jo N ^ nil’'*! ' ■ *•"’ 

^ t J L 


N N <*> 


-2 i I y v(f - f ) V (f - f - f ) 

N „-l n-l jio :l J 1 en 


N N 


4 il il 

* -2 I I V (f - f - f ) 

N 1 1 i 

“ m*l n“l 


N N 


= -2 I I V (f - (f - f )), 

N 0-1 n-l ^ 


which yields 


N-l N 


IPSD(£j)*IPSD(£^) = 1, I I V (fj - |f^„ - 

“ 0-1 n-o +1 


N-l N 


+ ^2 I I V (f + |f - f I) 

N 0-1 n-o +1 ^ 


+ ^aN V (f^) 


But > 0. Therefore, 


(D-10) IPSD(f.)*IPSD(f ) -^21 I 7 (f - 

m“l n*m+l 

Comparing Equation D-9 with D-10 gives 


(D-lla) IPSDg(f^) » irr-y IPSD(f^)*IPSD(f^) 


(D-llb) = I IPSD(f,)IPSD(f, - f,); 

N - 1 J J ^ 

i * all integers from 1 to + ®. 

Thus, given the IPSD, it is a relatively simple matter to 
compute the IPSDg. 

We are interested in the Averrgc Power Spectral Distribution 
of the beat frequency signals (APSD^) . But the average to be 
determined is the statistical average and not the time average. 

■y *T* 

(It can be shown that ^ f IPSD(fj^)dt' does not necessarily 
equal APSD_{f.).) Specifically, then, we want the statistical 
expectation of IPSDg(f^) at any time of day, t', with respect 
to the time of day, tj^2» etc., that the satellites can 

be at their perigees. This is mathematically expressed as 

(D-12) APSDg(f^) « E^j{lPSDg(f^) }. 

Using Equation D-il, Equation 0"12 becomes: 


APSD^(f^) - E^j 



I lPSD(f JIPSD(f 
j-0 - ^ 



> 0 . 


Using the expression for IPSD from Equation C-6 yields 


APSD 
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from Equation 10, p(tj^^) = 1/T, 0 < < T, all n. 
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The last term in Equation D-13 is the value of APSDg(fj^) when 

m « n in the summation over m and n. In this term, the 

product 6(f* - f ) 6(£" - f ) can only have a non-sero 

en en 

value when f ■ f « range of f is f^ ^ f* 

< fj + Lf} that of f" is fj - f^ f" £ f j + Af - fj^, 

and the two ranges do not have any values in common unless 
f^ = 0. However, ^ Equation D-13. Thus the last 

term in this equation is zero and we can drop the restriction 

m 7 ^ n on the other terms because m « n produces a zero value. 

Using the results of Equations C-13, C-17b, C-17c, C-18, 
and C-19, APSD_(f. ) in Equation D-13 becomes 

D 1 


APSDg(f j) 
(D-14) 


^ I APSD(f,)APSD(f 
N-1 j.o J 


* ^ APSD(f^)*APSD(f^) 






V fi > 0 




Thus, given the APSD(f^), it is a relatively simple matter to 
compute the APSDg via self-correlation of the APSD. (It is 
easy to show that the additional relationships below also hold 


Syff) 


2K*S (f) • S f) 

2K S(f')S(f 






f > 0 


f)df'J 


- 92 - 


Sb<‘> • 


f OO ^ A 

r\ 


2 K 


S(f')S(f - f)df 


y f > 0 


APSDg(f^) 


N(N-l) (KP)' 


+ Af 


Sg(f)df 


- 93 " 





